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Abstract 

Observations of the multi-TeV spectra of the nearby BL objects Mkn 421 and Mkn 
501 exhibit the high energy cutoffs predicted to be the result of intergalactic annihi- 
lation interactions, primarily with infrared photons having a flux level as determined 
by various astronomical observations. After correction for this absorption effect, 
the derived intrinsic spectra of these multi-TeV sources can be explained within 
the framework of simple synchrotron self-Compton emission models. Stecker and 
Glashow have shown that the existence of such annihilations via electron-positron 
pair production interactions up to an energy of 20 TeV puts strong constraints 
on Lorentz invariance violation. Such constraints have important implications for 
quantum gravity models as well as models involving large extra dimensions. We also 
discuss the implications of observations of high energy 7-rays from the Crab Nebula 
on constraining quantum gravity models. 
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1 Introduction 



In this paper, I will first present and discuss the evidence for absorption fea- 
tures at the multi-TeV end of the 7-ray spectra of the two active galactic 
nuclei for which we have the most detailed and highest quality spectral data. 
These 7-ray sources are the BL Lac objects known as Markarian (Mkn) 501 
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and Mkn 421. I will then show how this absorption is due to the expected 
7-ray annihilation caused by electron-positron pair production interactions of 
these 7-rays with intergalactic low energy photons having the flux level as 
determined by theoretical considerations combined with various astronomical 
observations. 

Stecker and Glashow have shown that the very existence of these interactions 
puts quantitative constraints on Lorentz invariance violation at high energies 
[1]. I will show here that this also implies significant constraints on a class 
of propsed models where quantum effects at the Planck length scale of ~ 
10~ 35 m, or a much larger modified Planck length scale in the case of large 
extra dimension models, can alter the relativistic energy-momentum dispersion 
relations for particles at energy scales much lower than the Planck scale (~ 
10 19 GeV). I will also discuss more severe indirect constraints on quantum 
gravity models which follow from observations of the 7-ray emission spectrum 
of the Crab Nebula. 



2 Absorption of Gamma-Rays at Low Redshifts 



The formulae relevant to absorption calculations involving pair-production are 
given and discussed in Ref. [2]. For 7-rays in the TeV energy range, the pair- 
production cross section is maximized when the soft photon energy is in the 
infrared range. In terms of the observed soft photon wavelength and 7-ray 
energy as a function of redshift, z, 



where A e = h/(m e c) is the Compton wavelength of the electron. For a 1 
TeV 7-ray emitted by a source at low redshifts, this corresponds to a soft 
photon having a wavelength in the J-band of IR astronomy (1.25/im). (Pair- 
production interactions actually take place with photons over a range of wave- 
lengths around the optimal value as determined by the energy dependence of 
the cross section.) If the emission spectrum of an extragalactic source extends 
beyond 20 TeV, then the extragalactic infrared field should produce an ab- 
sorption effect in the observed spectrum between ~ 20 GeV and ~ 5 TeV, 
depending on the redshift of the source [3]. 

Absorption of high energy 7-rays from extragalactic sources occurs via inter- 
actions of these photons with low energy photons of intergalactic radiation. 
These low energy photons are produced by stellar radiation and the reemission 
of such radiation by interstellar dust in galaxies. These photons then leave the 
galaxies in which they were produced, escaping into intergalactic space. 
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Fig. 1. The spectral energy distribution (SED) of the extragalactic infrared back- 
ground^]. The data are given along with the two semi- empirical model SEDs. The 
point at 100/im indicated by the small square derived in Ref. [11] is less certain than 
the data at lower frequencies because isotropy was not established for this point. 

Malkan and Stecker (MS01) [4] used empirically based spectral energy distri- 
butions (SEDs) of galaxies as a function of galaxy luminosity together with 
galaxy luminosity distribution functions to derive the SED of the cosmic in- 
frared background (CIB). The advantage of using empirical data to construct 
the SED of the CIB, as done in MS01, is particularly important in the mid IR 
range. In this region of the spectrum, galaxy observations indicate more flux 
from warm dust in galaxies than that taken account of in some more theoreti- 
cally oriented models, e.g., Ref. [5]. De Jager and Stecker (DS02) [6] extended 
these SEDs into the optical and UV using a hybrid model based on Hubble 
Space Telescope (HST) galaxy counts from the Hubble Deep Field (HDF). 
The lower curve in Figure 1 from DS02 (adapted from MS01) assumes a red- 
shift (temporal) luminosity evolution of galaxies oc (1 + z) 3 out to Zfi at = 2 
(baseline model), whereas the upper curve assumes "fast" oc (1 + z) A evolu- 
tion out to Zfi a t = 1.3. For z > Zfi at no evolution is assumed to occur. This is 
because evolution in star formation and stellar emissivity in galaxies appears 
to level off at redshifts greater than 2. [7]-[10] 

Figure 1 shows the SED curves from SD02 in comparison with various data and 
limits. The results of the MS01 [4] fast evolution SED generally agree well with 
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directly measured COBE (Cosmic Background Explorer) data. These results 
are also in agreement with upper limits obtained from TeV 7-ray studies [12] 
- [15] and lower limits obtained from various infrared galaxy count studies. 



3 Observations of the Blazars Mkn 501 and Mkn 421 

The highest energy extragalactic 7-ray sources in the known universe are the 
active galaxies called 'blazars', objects that emit jets of relativistic plasma 
aimed directly at us. Those blazars known as X-ray selected BL Lac objects 
(XBLs), or alternatively as high frequency BL Lac objects (HBLs), are ex- 
pected to emit photons in the multi-TeV energy range, but only the nearest 
ones are expected to be observable at TeV energies, the others being hidden 
by intergalactic absorption [16]. 

Extragalactic photons with the highest energies yet observed originated in a 
powerful flare coming from the giant elliptical active galaxy known as Mkn 501 
[17]. Its spectrum is most easily understood and interpreted as manifesting the 
high energy absorption to be expected from 7-ray annihilation by extragalac- 
tic pair production interactions. The analyses of de Jager and Stecker (DS02) 
[6] and Konopelko et al. [18] indicate the presence of the absorption effect 
predicted by calculating the expected energy dependent opacity inferred from 
the background light SEDs of of MS01 and DS02 (see previous section). This 
absorption is the result of electron-positron pair production by interactions of 
the multi-TeV 7-rays from Mkn 501 primarily with intergalactic infrared pho- 
tons. Intrinsic absorption by pair production interactions within 7-ray sources 
such as Mkn 421 and Mkn 501 is expected to be negligible because such giant 
elliptical galaxies contain little dust to emit infrared radiation and because 
BL Lac objects have little gas (and therefore most likely little dust) in their 
nuclear regions. It also appears that 7-ray emission in blazars takes place at 
superluminal knots in their jets, downstream of the radio cores of these active 
galaxies and therefore downstream of any putative accretion disks [19]. 

The spectrum of Mkn 501 in the flaring phase extends to an energy of at 
least 24 TeV [17]. The DS02 calculations predict that intergalactic absorption 
should strongly supress the spectra of these sources at multi-TeV energies. 
Figure 2 shows observed spectrum from HEGRA [17] and the Whipple tele- 
scope [20] (lower curve and points) and the derived intrinsic spectrum (upper 
curve and points) for Mkn501 in the flaring phase as given by de Jager and 
Stecker [6]. The intrinsic spectrum was derived by correcting for the opacity 
calculated for z = 0.03 as a function of energy, based on models of MS01, 
extended into the optical and UV range [6]. 

Konopelko, et al. [18] have reexmained the spectra of both Mkn 421 and Mkn 
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Fig. 2. The observed and derived intrinsic spectrum for Mkn 501 during the April 
1997 flare. The intrinsic spectra are given for the baseline and fast evolution models 
shown in Figure 1 [6]. 



501 corrected for absorption and have found that they fit synchrotron self- 
Compton (SSC) emission models from the radio to TeV 7-ray range. Both 
the X-ray and intrinsic TeV spectra of Mkn 421 peak at lower energies than 
those of the flaring spectrum of Mkn 501 (see Figs. 3 and 4) which can be 
understood if electrons were accelerated to higher energies in the April 1997 
flare of Mkn 501 than in Mkn 421. As shown in Figure 2, the intrinsic spectrum 
of Mkn 501 with the absorption effect removed actually peaks at multi-TeV 
energies, rather than falling off in this energy range. Thus, it appears that the 
dropoff in the observed 7-ray spectrum of Mkn 501 above ~ 5 TeV is a direct 
consequence of intergalactic absorption. We will therefore interpret the Mkn 
501 data as evidence for intergalactic absorption with no indication of Lorentz 
invariance violation (see next section) up to a photon energy of ~ 20 TeV. This 
conclusion can be substantiated by observations of other extragalactic 7-ray 
sources at differing redshifts [3]. 



4 High Energy Consequences of Breaking of Lorentz Invariance 



It has been suggested that Lorentz invariance (LI) may be only an approxi- 
mate symmetry of nature [21] - [23]. A simple and self-consistent framework 
for analyzing possible departures from exact LI was suggested by Coleman 
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Fig. 3. The April 1997 flaring spectrum of Mkn501 with absoprtion taken into 
account at TeV energies, shown along with the SSC model fit of Konopelko et al. 
[18] over a large energy range. 



and Glashow [24], who assume LI to be broken perturbatively in the con- 
text of conventional quantum field theory. Small Lorentz noninvariant terms 
are introduced that are renormalizable, being of mass dimension no greater 
than four and having dimensionless coupling constants in the kinetic part 
of the Lagrangian. These terms are also chosen to be gauge invariant under 
£77(3) x SU{2) x U(l) ("the almost standard model"). It is further assumed 
that the Lagrangian is rotationally invariant in a preferred frame which is 
presumed to be the rest frame of the cosmic microwave background. 

Consequent observable manifestations of LI breaking can be described quite 
simply in terms of different maximal attainable velocities of different particle 
species as measured in the preferred frame. This is because the small LI vio- 
lating terms modify the free-field propagators so that the maximum velocities 
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Fig. 4. The spectra of Mkn421 in two phases with absoprtion taken into account at 
TeV energies, shown along with the SSC model fit of Konopelko et al. [18] over a 
large energy range. 

of various particles are not equal to c. Indeed, this type of LI breaking within 
the hadron sector is one way to circumvent the predicted but unseen 'GZK 
cutoff' in the ultrahigh energy cosmic-ray spectrum owing to photomeson in- 
teractions of such cosmic rays with photons of the 2.7K cosmic background 
radiation [25], [26]. These interactions are expected to produce an effective 
attenuation mean-free-path for ultrahigh energy cosmic rays (E > 10 20 eV) in 
intergalactic space of < 100 Mpc [27] in the absence of LI breaking. 



5 The LI Breaking Parameter S 

The discussion in the previous section shows that if LI is violated the maximum 
attainable velocity of an electron need not equal the velocity of light in vacuo, 
i.e., c e 7^ c 7 . The physical consequences of this violation of LI depend on 
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the sign of the difference between c e and c 7 [24]. Following the discussion of 
Stecker and Glashow [1], we define 



and consider the two cases of positive and negative values of 5 separately. 

Case I: If c e < c 7 (5 < 0), the decay of a photon into an electron-positron pair 
is kinematically allowed for photons with energies exceeding 



The decay would take place rapidly so that photons with energies exceeding 
E max could not be observed either in the laboratory or as cosmic rays. Since 
photons have been observed with energies E 1 > 50 TeV from the Crab nebula 
[28], we deduce for this case that £ max > 50 TeV, or that \5\ < 2 x 10~ 16 
[1]. Stronger bounds on 5 can be set through observations of very high energy 
(TeV) photons. The detection of cosmic 7-rays with energies greater that 
50 TeV from sources within our galaxy would improve the bound on 5. 

Case II: Here we are concerned with the remaining possibility, where c e > c 7 
(5 > 0) and electrons become superluminal if their energies exceed E mSLX /\/2. 
Electrons traveling faster than light will emit light at all frequencies by a 
process of 'vacuum Cerenkov radiation.' This process occurs rapidly, so that 
superluminal electron energies quickly approach E max /^/2. Because electrons 
have been seen in the cosmic radiation with energies up to ~ 2 TeV[29], it 
follows that 5 < 3 x 10~ 14 . This upper limit is about two orders of magnitude 
weaker than the limit obtained for Case I. 

A smaller, but more indirect, upper limit on 5 for the 5 > case can be ob- 
tained from theoretical considerations of 7-ray emission from the Crab Nebula. 
The unpulsed 7-ray spectrum of the Crab Nebula can be understood to be 
produced by sychrotron emission up to the ~ 0.1 GeV 7-ray energy range (see 
Figure 5). Above 25 MeV, the synchrotron component falls off very rapidly 
with energy as expected from theoretical limits on electron acceleration [30]. 
Emission above 0.1 GeV, extending into the TeV range, can be explained as 
synchrotron self-Compton emission of the same relativistic electrons which 
produce the synchrotron radiation. The Compton component, extending to 
50 TeV, implies the existence of electrons having energies at least this great 
in order to produce 50 TeV photons, even in the extreme Klein-Nishina limit. 
This is, of course, required by conservation of energy. This indirect argument, 
based on the reasonable assumption that the 50 TeV 7-rays are from Comp- 
ton interactions, leads to a smaller upper limit on 5, viz., 5 < 10~ 16 . Better 
observational 7-ray data for the Crab Nebula at multi-TeV energies, and data 



c e = c 7 (l + 5) , < \S\ < 1 , 



(2) 
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Fig. 5. The observed 7-ray spectrum of the Crab nebula with curves showing the 
synchrotron and Compton components (from Ref. [30]). 

on the shape of its 7-ray spectrum above 50 TeV, will be needed to further 
test this indirect constraint on 5. One open question involves the possibility 
of hadronic interactions involving cosmic ray nucleons producing 7T°'s which 
decay into observed high energy 7-rays(although there is no present indica- 
tion of this). Hopefully, better data will determine whether or not there is 
a significant hadronically induced component contributing to the multi-TeV 
spectrum of the Crab. 

A further constraint on 5 for 5 > (c e > c 7 ) follows from the modification of 
the threshold energy for the pair production process 7 + 7 — > e + + e~. This 
follows from the fact that the square of the four-momentum is changed to give 
the threshold condition 



where e is the energy of the low energy photon and 6 is the angle between the 
two photons. The second term on the left-hand-side comes from the fact that 
c 7 = dE 1 /dp 1 . 



2e£ 7 (l - cos 9) - 2E 2 5 > 4m 2 e , 



(4) 
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For head-on collisions (cos# = —1) the minimum low energy photon energy 
for pair production becomes 




m 2 e E^5 



(5) 



mm 



2 



It follows that the condition for a significant increase in the energy threshold 
for pair production is E^S/2 > m^/Ey, or equivalently, 



As discussed in the previous section, there is no indication of LI violation 
supressing the physics of pair production for photons up to an energy of ~ 20 
TeV. Thus, it follows from eq. (6) that the Mkn 501 observations imply the 
constraint 5 < 2ml/ E 2 = 1.3 x 10~ 15 [1]. This constraint on positive 5 is 
more secure than the smaller, but indirect, limit given by the Crab Nebula 
acceleration model. 



6 Quantum Gravity Models 

In the absence of a true and complete theory of quantum gravity, theorists 
have been suggesting and exploring models to provide experimental and ob- 
servational tests of possible manifestations of quantum gravity phenomena. 
Such phenomena have usually been suggested to be a possible result of quan- 
tum fluctuations on the Planck scale M P i anck = y/he/G ~ 1.22 x 10 19 GeV/c 2 , 
corresponding to a length scale ~ 1.6 x 10~ 35 m [31] - [33]. In models involv- 
ing large extra dimensions, the energy scale at which gravity becomes strong 
can be much smaller than Mpi anc k, with the quantum gravity scale, Mqg, 
approaching the TeV scale [34], [35]. 

In many of these models Lorentz invariance is predicted to be violated at 
high energy. This results in interesting modifications of particle physics that 
are accesible to observational tests using TeV 7-ray telescopes and cosmic ray 
detectors. An example of such a model is the loop quantum gravity model 
with a preferred inertial frame given by the cosmological rest frame of the 
cosmic microwave background radiation (For an extensive discussion, see the 
review given in Ref. [36].) 

In the most commonly considered of these models, the usual relativistic disper- 



S > 




(6) 



10 



sion relations between energy and momentum of the photon and the electron 



E> = P* (7) 
El = Pl + m 2 e (8) 

(with the "low energy" speed of light, c = 1) are modified by a leading order 
quantum space-time geometry corrections which are cubic in p ~ E and are 
supressed by the quantum gravity mass scale Mqg- Following Refs. [22] and 
, we take the modified dispersion relations to be of the form 

E 2 = v 2 ^— (9) 



El + (10) 

M QG 
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We assume that the cubic terms are the same for the photon and electron 
as in eqs. (9) and (10). More general formulations have been considered by 
Jacobson, Liberati and Mattingly [37] and Konopka and Major [38]. 

As opposed to the Coleman- Glashow formalism, which involves mass dimen- 
sion four operators in the Lagrangian and preserves power-counting renor- 
malizablility, the cubic term which modifies the dispersion relations may be 
considered in the context of an effective "low energy" field theory, valid for 
E <C Mqg, in which case the cubic term is a small perturbation involving di- 
mension five operators whose construction is discussed in Ref. [39]. With this 
caveat, we can generalize the LI violation parameter 5 to an energy dependent 
form 



dE e dE 1 ^ E 7 ml E e 

~dp~ ~" ~dp^ ~ Mq^ ~ 2Ej " Mq? ( } 



which is a valid approximation for the energy regime E e ^> m e . Note that the 
maximum velocities of particles of type % are reduced by 0{Eij 'Mqg). 

For pair production then, with the positron and electron energy E e ~ E y /2, 

X _ ^ 2 ™e (12) 



2M QG E* 



and the threshold condition given by eq.(6) reduces to the constraint 

M QG > g. (13) 
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Since pair production occurs for energies of at least 20 TeV, as indicated 
by our analyses of the Mkn 501 and Mkn 421 spectra [6], [18], we then find 
the constraint on the quantum gravity scale Mq G > 0.3M Planck . This cons- 
triant contradicts the predictions of some proposed quantum gravity models 
involving large extra dimensions and smaller effective Planck masses. Pre- 
vious constraints on Mqq for the cubic model, obtained from limits on the 
energy dependent velocity dispersion of 7-rays for a TeV flare in Mkn 421 
[40] and from 7-ray bursts [41] were in the much less restrictive range Mqq > 
(5 - 7) x 10^M Planck . 



7 Beyond the Planck Scale: Implied Constraints from Crab Nebula 
7-rays 

Within the context of a more general cubic modification of the dispersion 
relations given by eqs. (9) and (10), Jacobson, et al. [42] have obtained an 
indirect limit on Mqq from the apparent cutoff in the synchrotron component 
of the in the Crab Nebula 7-ray emission at ~ 0.1 GeV (see Figure 5). By 
making reasonable assumptions to modify the standard synchrotron radiation 
formula to allow for Lorentz invariance violation, they have concluded that 
the maximium synchrotron photon energy will be given by -E 7)max = 0.34 
(eB /m e )(m e / Mqg)~ 2 ^ ■ This reasoning leads to the constraint Mq G >1.2 x 

10 M Planck- 

Future observations of the Crab Nebula with the GLAST (Gamma-Ray Large 
Area Space Telescope) satellite, scheduled to be launched in 2005, will provide 
a better determination of its unpulsed 7-ray spectrum in the energy range 
above 30 MeV where the transition from the synchrotron emission component 
to the Compton emission component occurs. This will provide a more precise 
determination of the maximum electron energy in the Nebula and therefore 
provide a more precise constraint on the parameter Mq G as we have defined 
it here. However, this constraint will still be orders of magnitude above the 
Planck scale. 



8 Conclusions 

Nearly a century after the inception of special relativity, high energy 7-ray 
observations have confirmed its validity up to electron energies of 2 TeV, 
photon energies of 20 TeV and, indirectly, up to electron energies in the PeV 
range. These results indicate an absence of evidence for proposed violations of 
Lorentz invariance as predicted by some phenomenological quantum gravity 
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and large extra dimension models. Thus, high energy astrophysics has provided 
important empirical constraints on Planck scale physics. 

Models with large extra dimensions are ruled out by the existence of absorption 
in the very high energy spectra of nearby BL Lac objects. The fact that more 
distant brighter sources are not seen can also be taken as indirect evidence of 
intergalactic absorption by pair production interactions [16]. 

The constraints based on analysis of the Crab Nebula 7-ray spectrum, dis- 
cussed in the previous section, imply that the quantum gravity scale is orders 
of magnitude above the Planck mass scale. This indicates that the class of 
models considered here cannot be reflective of physics at the Planck scale. 
Models such as loop quantum gravity with a preferred inertial frame are ruled 
out by this line of reasoning. Alternative models to consider might be models 
with a quartic term with Mq G supression in the dispersion relations, Lorentz 
invariant quantum gravity models, or really new Planck scale physics such as 
string theory, which preserves Lorentz invariance. 
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